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is  sometimes  able  to  discriminate  bodies  having  high  values  of  form  drag. 

It  is  concluded  that  the  simple  drag  formula  may  not  be  reliably  used  for 
estimating  the  relative  form  drags  of  bodies  of  revolution. 
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ABSTRACT 

A simplified  drag  formula  used  by  Hess  and  James  to  predict  the 
drag  of  a small  number  of  axisymmetric  forms  is  systematically  applied 
to  seven  series  of  model  forms  comprising  nearly  fifty  bodies. 

Calculations  of  form  (or  residual)  drag  are  compared  to  available 
experimental  data  in  order  to  determine  the  usefulness  of  the  method 
for  predictive  purposes.  The  formula  is  shown  to  exhibit  very  little 
sensitivity  to  changes  in  most  body  parameters  such  as  the  length-of- 
stern  to  diameter  ratio,  nose  and  tail  radii  parameters,  and  the  prismatic 
coefficient.  For  some  parameters,  such  as  length  to  diameter  ratio, 
and  length  of  bow  section  to  diameter  ratio,  however,  the  formula  is 
sometimes  able  to  discriminate  bodies  having  high  values  of  form  drag. 

It  is  concluded  that  the  simple  drag  formula  may  not  be  reliably  used 
for  estimating  the  relative  form  drags  of  bodies  of  revolution. 

ADMINISTRATIVE  INFORMATION 

This  work  was  authorized  and  funded  by  the  Naval  Sea  Systems  Command, 

Task  Area >^66007,  Task  08101,  Work  Unit  1552-135. 


INTRODUCTION 
1 2 

In  a recent  report  and  paper,  Hess  and  James  compare  the  predictions 
of  a simple  drag  formula  with  experimental  data  for  eight  Series  58  bodies 
of  revolution.  The  small  number  of  bodies  used  by  Hess  and  James  to  justify 
their  use  of  this  simple  formula  is  deemed  insufficient  to  establish  the 
method's  validity.  The  present  investigation  attempts  to  remedy  this  by 
making  systematic  tests  of  the  formula  on  many  more  axisymmetric  bodies 
for  which  drag  data  are  available.  This  study  attempts  to  determine  whether  or 


AHess,  John  L.  and  R.N.  James,  "On  the  Problem  of  Shaping  an  Axisymmetric 
Body  to  Obtain  Low  Drag  at  Large  Reynolds  Numbers,"  McDonnell  Douglas 
-Corporation  Report  MDC-J6791  (Jan  1975)  . 

^Hess,  John  L.,  "On  the  Problem  of  Shaping  an  Axisymmetric  Body  to  Obtain 
Low  Drag  at  Large  Reynolds  Numbers,"  Journal  of  Ship  Research,  Vol.  20, 
No.  1 (May  1976). 

uertler,  Morton,  "Resistance  Experiments  on  a Systematic  Series  of 
Streamlined  Bodies  of  Revolution  - For  Application  to  the  Design  of  High- 
Speed  Submarines,"  David  Taylor  Model  Basin  Report  C-297  (Apr  1950), 
declassified  27  Jan  1967. 
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not  the  simple  formula  can  predict  the  same  order  of  "drag"  merit  as  that 
measured  experimentally  for  a series  of  bodies  of  revolution.  Even  if  the 
drag  coefficients  are  not  in  excellent  numerical  agreement,  the  formula 
might  be  useful  as  a relative  indicator  of  drag  characteristics. 

THE  SIMPLE  DRAG  FORMULA 

The  drag  of  a body  is  usually  calculated  in  one  of  two  ways:  (1)  by 

integration  of  local  surface-pressure  and  skin-friction  forces  acting  on  the 

body  surface,  or  (2)  by  determination  of  the  velocity  profile  and  momentum 

defect  in  the  wake  far  downstream  of  the  body  where  the  pressure  is  equal  to 

4 5 

the  ambient  pressure.  Young  and  Granville  provide  outlines  of  the 
derivation  of  the  momentum  integral  for  the  latter  method  which  gives  the 

total  drag,  D,  of  a body  of  revolution  as 

00 

D =*  2irp  J u (U^  - u)  r dr  [l] 

0 


where  ■ velocity  of  incoming  flow 

u = velocity  of  the  fluid  in  the  wake 
r = radial  distance  from  the  x axis 
p = fluid  density 

The  general  body  orientation  and  parameter  notation  can  be  found  in 
Figure  1. 

If  we  define  the  momentum  area  ft  as 

[2] 

where  6 is  the  half-thickness  of  the  wake,  then,  for  the  wake  far 
downstream  ft  = ft  , with  <5  +*,  and 

4Young,  A.D.,  "The  Calculation  of  the  Total  and  Skin  Friction  Drag  of 
Bodies  of  Revolution  at  Zero  Incidence,"  Aeronautical  Research  Committee 
(.R6M  Report  1874  (Apr  1939)  . 

^Granville,  P.S.,  "The  Calculation  of  the  Viscous  Drag  of  Bodies  of 
Revolution,"  DTMB  Report  849  (Jul  1953). 
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Figure  1 - Body  Orientation 
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Hence,  far  a drag  coefficient,  C , based  on  wetted  surface  area  A 
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[4] 
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Young  indicated  that  the  wake  downstream  of  a body  is  actually  a 
continuation  of  the  boundary  layer  on  the  body  tail*  Thus,  boundary  layer 
theory  may  be  used  to  determine  the  momentum  area,  ft,  at  the  end  of  the 
tail.  Then,  it  is  necessary  to  relate  ft  at  the  body  tail  to  ft^  far  down- 
stream. 

For  turbulent  flow,  the  equations  of  motion  are  given  by  the  Reynolds 
equations  (the  averaged  Navier-Stokes  equations) . For  axisymmetric  flow 
past  a body  of  revolution  having  sufficiently  small  longitudinal  curvature, 
the  equations  of  motion  from  Reference  5 are 
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where  x,  y,  $ are  body  of  revolution  coordinates  (see  Figure  1)  with 
average  velocity  components  u,  v,  and  w = 0 and  turbulent  velocity 
components  u',  v'  and  w'. 


The  continuity  equation  is  given  by 


3(ru)  3(rv) 

+ 0 [8] 

3x  3y 

Now,  using  Prandtl's  assumption  that  the  boundary  layer  is  thin,  and 

2 3 

assuming  that  the  Reynolds  number  is  of  order  1/6  , 3/3x  <<  3/3y,  — 'v  1/6 
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8/ 


and  — 2 % such  that  v'  is  of  order  6 and  u,  u',  3/3x  and  r are  of  order 

0/  6Z 

1,  equations  [5-7]  reduce  to 
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with 
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3y 


For  additional  information  on  the  process  of  reduction  see  Reference  6. 
Integration  of  equation  [9]  across  the  boundary  layer  from  y = 0 to  y = 6, 
with  use  of  the  continuity  equation  [8]  and  the  definition  of  equation  [2], 
gives  the  momentum  equation  for  axisymmetric  turbulent  boundary  layers  on 
bodies  of  revolution: 


ft 


1 Jrdy 

U J 


[10] 


where  h =- 


the  axisymmetric  shape  factor 


j “ 

n u. 
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p pu'v' , shear  stress 

3y 

a = p (u')2 


r , t = values  of  r and  t at  the  wall, 
w w 

U = velocity  at  the  edge  of  the  wake 


Schlichting,  H.,  Boundary  Layer  Theory,  McGraw  Hill,  New  York,  1968. 
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If  it  is  assumed  that  (u‘)  /u<<1,  then  for  wake  flow  having  no  wall  shear 
this  reduces  nicely  to 


3n  n au 

— + (h  + 2)  — — - 0 

3x  U 3x 


The  conditions  at  the  end  of  the  body  (denoted  by  subscript  "e") 
can  be  related  to  conditions  downstream  in  the  wake  by  integrating  [11] 
by  parts.  This  yields 


n - n 
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I ) exp  I In  — dh 
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00  LI 


Several  experimental  curve  fittings  for  this  equation  have  been  made, 
among  them  fits  by  Young  and  Granville.  Granville  suggest  using 
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where  q is  determined  experimentally.  If  this  is  substituted  into  [12], 


there  results 
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which  when  substituted  into  [4]  yields 


(he+2)q+3 


C_  = n 

D . e ' 
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To  obtain  12  it  is  necessary  to  return  to  equation  [10]  . Under  the 

e 7 

same  assumptions  made  in  deriving  equation  [11],  Granville  develops  a 

relation  for  bodies  of  revolution  having  thin  boundary  layers, 

7Granville,  P.S.,  "Partial  Form  Factors  from  Equivalent  Bodies  of 
Revolution  for  the  Froude  Method  of  Predicting  Ship  Resistance,"  SNAME 
First  Ship  Technology  and  Research  Symposium,  Washington,  D.C.  (Aug  1975). 
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where  the  derivative  is  in  terms  of  ds  instead  of  dx  (dx  = sec  a ds) . Here 


s = distance  along  the  body  meridian 

r = the  radial  distance  of  the  body 
w 

t*  = wall  shear  stress  = t sec  a 
w w 

U = the  velocity  at  the  edge  of  the  boundary  layer. 


Granville  extends  the  formulation  to  thick  axisymmetric  boundary  layers 
writing 


PU 


' w' 


where  k and  m are  constants,  and  v is  the  kinematic  viscosity  of  the  fluid. 
Then  assuming  that  h = constant  = hg  the  momentum  equation  [16]  may  be 
integrated  to  give: 


(l-hn)(h+2)  . v ( 1+m)  ( 1 >(14i)/d  . [(Hm)  (h+2)-m] 

) (?)  U —fe 


where 

L = body  length 

Rn.  = Reynolds  number  based  on  Length  = 

L v 

a = angle  between  the  body  axis  and  a tangent  to  the  meridian 
at  given  x 


Substitution  of  the  value  of  determined  from  [17]  into  [15]  gives 
the  drag  formula 
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For  a flat  plate  U = U^,  sec  a = 1,  and  rw  = A/2„l,  so  that  the 
drag  coefficient  CD  reduces  to  a frictional  drag  coefficient,  C^,  given 
bX 


[18] 
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which,  if  substituted  into  equation  [18]  gives 
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[19] 


where  Cf  is  the  frictional  drag  coefficient  of  a flat  plate  having  the 

same  length  and  Reynolds  number  as  the  body  of  revolution  being  considered. 

Now,  equation  [19]  is  a general  form  of  the  simple  drag  formula  derived  in 

7 2 

different  ways  by  Granville  and  Hess.  Values  of  m,  q,  and  h employed  by 
the  two  authors  are 


Granville:  m = 0.1686,  q = 7,  h = 1.4 
Hess:  m = 1/6,  q = 1,  h = 1 

It  should  be  noted  that  Hess's  numbers  were  derived  from  work  by 
8 9 

Truckenbrodt  and  Falkner.  For  this  report,  the  values  used  were  those 
chosen  by  Hess  for  his  work.  Thus,  formula  [19]  reduces  to 


[feTe 
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[20] 


g 

Truckenbrodt,  E.,  "A  Method  of  Quadrature  for  Calculation  of  the  Laminar 
and  Turbulent  Boundary  Layer  in  Case  of  Plane  and  Rotationally  Symmetrical 
qFlow,"  NACA  TM  Report  1379  (May  1955). 

'Falkner,  V.N.,  Aircraft  Engineering,  15  (1943),  p.  65. 
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OUTLINE  OF  COMPUTATIONS 


Point-generating  computer  programs  were  developed  to  give  an  accurate 
physical  representation  of  each  body  surface  for  computer  use.  The  creation 
of  these  point  distributions  usually  involved  using  the  original  model 
drawings  or  offsets  to  obtain  a fairly  accurate  curve.  A least-squares 
polynomial  fit  was  then  passed  through  the  data  points  to  insure  smooth  rates 
of  change  in  r.  Two  series  of  bodies  were  an  exception  to  this;  they  were 
generated  analytically  by  using  the  appropriate  functions  from  which  they 
were  developed.  The  points  generated  were  fed  into  an  existing  Douglas- 
Neumann  axisymmetric-body  program,  DASO,^  to  compute  the  distributions  of 
potential-flow  velocity  U/U,  cos  a,  and  pressure  coefficient,  C^. 

With  the  output  of  the  axisymmetric  program,  calculations  of  the 
integrand  (U/U^)^0^3  (r/L)^  sec  a,  appearing  in  formula  [20],  were  made 
and  graphed  against  x/L.  Typical  distributions  of  the  integrand  are  shown  in 
Figure  2.  The  solid  curve  illustrates  the  integrand  for  a model  with  parallel 
middle  body  and  the  broken  line  indicates  the  integrand  for  a typical  model 
having  no  parallel  middle  body.  Using  data  points  spaced  every  0.5  percent 
for  the  first  20  percent  of  the  body  and  every  1 percent  thereafter,  a fine 
grid  integration  was  performed  to  compute  C^/C^.  Spacings  this  fine  were 
required  to  achieve  accurate  computer  results. 

Since  most  publications  report  the  residuary-drag  coefficient  CR  which 
is  defined  by  CR  = C^  - C^.,  instead  of  C^,  CR/Cj-  was  also  calculated. 
Comparisons  between  computational  results  and  experimental  data  were  made 
on  the  basis  of  CR/C^  and  C.R. 

No  attempt  is  made  in  this  report  to  evaluate  the  accuracy  of  the 
available  experimental  data,  some  of  which  may  have  been  subject  to  error. 
Comparisons  are  made  only  on  relative  order  of  merit  as  found  by  the  various 
experimenters  and  as  computed  by  the  simplified  formula. 

RESULTS 

Seven  series  of  bodies  were  used  to  test  the  simple  drag  formula: 

(1)  the  Series  58  bodies  investigated  by  Gertler,3  (2)  Series  5242  sterns 

10Hess,  J.L.  and  A.M.O.  Smith,  "Calculation  of  Potential  Flow  About 
Arbitrary  Bodies,"  in  "Progress  in  Aeronautical  Sciences,"  Tergamon 
Press,  Oxford  and  New  York,  Vol . 8 (1966),  pp.  1-138. 
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Figure  2 - Typical  Curves  for 


(unpublished),  (3)  six  bodies  with  inflected  sinusoidal  sterns  investigated 

11  12 

by  Kempf  (herein  after  referred  to  as  "Kempf  bodies"),  * (4)  a series  of 

bodies  based  on  the  best  Series  58  model  with  various  amounts  of  parallel 

middle  body  as  reported  by  Larsen,15  (5)  Series  4620  forebodies  investigated 

14 

by  McCarthy,  Power,  and  Huang,  (6)  polynomial  representations  of  five 
miscellaneous  models  and  finally,  (7)  Series  4935  afterbodies  (unpublished). 


Series  58 

The  first  series  of  bodies  was  Series  58,  originally  developed  by 
Gertler5  and  extensively  studied  by  Landweber  and  Macagno.  This  series  is 
a family  of  bodies  of  revolution  with  sectional -area  curves  represented  by 
sixth-degree  polynomials.  The  coefficients  of  the  equations  depend  on  six 
geometric  quantities:  body  length,  L,  maximum  diameter,  d,  radii  of  curvature 

at  the  nose  and  tail,  R0  and  R^ , body  volume,  V,  and  axial  position  of  body 
diameter,  x^.  From  these  six  parameters,  several  nondimens ional  quantities 
were  formed  and  used  by  Landweber  and  Macagno  to  study  variations  in  body 
shape.  The  nondimensional  quantities  used  were  length-to-diameter  ratio, 

X = L/d,  position  of  maximum  section,  m = x /L,  nondimensional  nose  and  tail 

2 2 m 2 
radii,  r0  = RQL/d  and  = R2L/d  , and  the  prismatic  coefficient,  = 4V/nLd  . 

Nondimensional  potential  flow  velocities  on  the  surfaces  of  the  Series  58 

bodies  were  published  by  Landweber  and  Macagno15  and  were  used  here  in  computing 

the  values  of  C^/C^.  by  the  simple  formula,  equation  [20].  Sketches  of  the 

bodies  can  be  found  in  reference  15  as  can  the  computed  pressure  distributions. 

11 Kempf,  George,  "Resistance  and  Wake  of  Some  Bodies  of  Revolution,"  from 
"New  Developments  in  Ship  Research,"  Jahrbuch  Schiffbautechnischen 
j^Gesel lschaft  (1927),  pp.  177-178. 

Kempf,  George,  "Turbulent  Separation  on  Full  Ship  Forms,"  Schiff  und  Hafen, 
j-Vol . 6,  No.  7,  Hamburg  (1954). 

Larsen,  C.A.,  "Additional  Tests  of  Series  58  Forms,  Part  1,  Resistance  Tests 
of  a Parallel  Middle  Body  Series,"  David  Taylor  Model  Basin  Report  C-738 
j^(Nov  1955),  declassified  on  2 Sep  1975. 

1 McCarthy,  J.H.,  J.  Power,  and  T.T.  Huang,  "The  Roles  of  Transition,  Laminar 
Separation  and  Turbulence  Stimulation  in  the  Analysis  of  Axisymmetric  Body 
Drag,"  to  be  published  in  the  Proceedings  of  the  Eleventh  ONR  Symposium  on 
Naval  Hydrodynamics,  sponsored  by  the  Office  of  Naval  Research,  London, 

(Mar  1976). 

15Landweber,  L.  and  Matilde  Macagno,  "Potential  Flow  about  Series  58  Bodies 
in  General  Translational  and  Rotational  Motion,"  Naval  Ship  Research  and 
Development  Center  Report  2505  (Jun  1967) . 
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TABLE  1 


I 

r 


CD  Comparisons  for  Series  58 

K 


Body 

No. 

Experimental 
CR  x 103 

Calculated 

Vcf 

Simple 

Formula 

Calculated 

C * 3 

lR  x 10 

Experimental 
Order  of  merit 

Computed 
order 
of  merit 
simple 
formula 

4154 

0.58 

0.184 

0.48 

24 

24 

4155 

0.36 

0.141 

0.37 

21 

23 

4156 

0.22 

0.115 

0.30 

18 

20 

4157 

0.13 

0.089 

0.23 

6 

16 

4158 

0.09 

0.068 

0.18 

3 

3 

4159 

0.075 

0.043 

0.11 

2 

1 

4160 

0.12 

0.078 

0.21 

5 

5 

4161 

0.15 

0.080 

0.21 

12 

6 

4162 

0.17 

0.084 

0.22 

15 

8 

4163 

0.19 

0.084 

0.22 

17 

9 

4164 

0.37 

0.092 

0.24 

22 

17 

4165 

0.07 

0.088 

0.23 

1 

14 

4166 

0.28 

0.092 

0.24 

19 

18 

4167 

0.16 

0.094 

0.25 

13 

19 

4168 

0.14 

0.085 

0.22 

10 

11 

4169 

0.14 

0.088 

0.23 

11 

15 

4170 

0.18 

0.077 

0.20 

16 

4 

4171 

0.13 

0.063 

0.17 

7 

2 

4172 

0.13 

0.082 

0.22 

8 

7 

4173 

0.13 

0.085 

0.22 

9 

12 

4174 

0.10 

0.084 

0.22 

4 

10 

4175 

0.32 

0.137 

0.36 

20 

22 

4176 

0.41 

0.133 

0.35 

23 

21 

4177 

0.16 

0.087 

0.23 

14 

13 

Based  on  Rn^=  20  x 10^ 

Cf  = 2.628  x 10’3  (Schoenherr1^) 
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For  easy  reference,  the  computed  pressure  distributions  are  included  here 
in  Appendix  A,  Figures  la  through  7A. 

Table  1 shows  the  results  for  the  Series  58  bodies  ordered  by  model 


number.  The  order  of  merit  shown  is  measured  relative  to  the  model  having 
the  smallest  experimental  value  of  CD.  It  is  apparent  that  for  the  overall 
series,  the  computed  orders  of  merit  are  quite  different  from  the 
experimentally  determined  orders  of  merit.  The  best  body  according  to 
Gertler's  experiments  was  4165.  This  is  No.  14  in  order  of  merit  according 
to  the  simple  formula.  The  only  agreement  is  in  the  relative  placement  of 
bodies  4158,  4160,  and  4154.  The  calculated  values  of  CR  shown  in  Table  1 
were  computed  from  CD/Cf  using  the  constant  value  of  Cc  appropriate  to  a 
Reynolds  number  of  20  x 10  ; this  value  of  C^.  was  also  used  by  Hess^  for 
his  comparisons. 

As  mentioned  above,  several  nondimensional  parameters  were  formed  from 
the  geometric  properties  of  the  Series  58  models.  These  nondimensional 
quantities  were  varied  systematically  to  create  several  subseries  within  the 
overall  model  range  covered.  Tables  2,  3,  4,  5,  and  6 list  the  bodies  involved 
in  these  seVies  and  compare  the  experimental  results  with  the  drag  predictions 
of  the  simple  formula.  Figure  5 provides  a graphical  comparison  between  the 
experimental  and  calculated  values  of  CD.  Points  are  labeled  in  terms  of  the 
major  parameter  varied,  and  not  by  body  number. 

TABLE  2 

Series  58  Body  Variations  in  L/D 


(m  = 

0.40,  C = 0.65, 
P 

rQ  = 0.5,  rx  = 0.1) 

Model 

Number 

L/D 

Experimental 
CR  x 103 

Gertler 

Computer 

S/Cf 

Simple  Formula 

Computed^ 

V x 10 

Simple  Formula 

4154 

4.0 

0.58 

0.184 

0.48 

4155 

5.0 

0.36 

0.141 

0.37 

4156 

6.0 

0.22 

0.115 

0.30 

4157 

7.0 

0.13 

0.089 

0.23 

4158 

8.0 

0.09 

0.068 

0.18 

4159 

10.8 

0.075 

0.043 

0.11 

*Based  on 

RnL=  20 

x 106  (16) 
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TABLE  3 

Series  58  Body  Variation  in  m,  Maximum  Thickness  Position 
(L/D  = 7.0,  C = 0.65,  ro  = 0.5,  ^ = 0.1) 


Model 

Experimental 
CR  x 103 

Computed 

Cn/Cf 

Computed 
c * 

R j 

Number 

m 

Gertler 

Simple  Formula 

Simple  Formula 

4160 

0.36 

0.12 

0.078 

0.21 

4161 

0.44 

0.15 

0.080 

0.21  { 

4162 

0.48 

0.17 

0.084 

0.22  j 

4163 

0.52 

0.19 

0.084 

0.22  j 

i 

} 

i 

| 

• 

TABLE 

4 

Series 

58  Body  Variations  in  C , 

The  Prismatic  Coefficient  ] 

(L/D 

= 7.0,  m = 0.4, 

rQ  = 0.5,  r1  » 0 

.i)  i 

Model 

r 

Experimental 
CR  x 103 

Computed 

CR/Cf 

Computed  j 

Cr* 

Number 

p 

Gertler 

Simple  Formula 

Simple  Formula 

4164 

0.55 

0.37 

0.092 

0.24 

4165 

0.60 

0.07 

0.088 

0.23 

4166 

0.70 

0.28 

0.092 

0.24 

"Based 

on  Rn^=  20 
Cf  = 2- 

x 106 

638  x 10‘3  (16) 
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TABLE  5 


Series  58  Body  Variations  in  r , Nondimens ional  Nose  Radius 
(L/D  = 7.0,  Cp  = 0.65,  m = 0.4,  rj  = 0.1) 

Experimental  Computed  Computed 

C * 


Model 

Number 


Experimental 
CR  x 103 

Gertler 


Vcf 


Simple  Formula  Simple  Formula 


0.094 

0.085 

0.088 

0.077 


TABLE  6 

Series  58  Body  Variations  in  r^ , Nondimens ional  Tail  Radius 
(L/D  = 7.0,  Cp  = 0.65,  m = 0.4,  rQ  = 0.5) 


Model 


Experimental  Computed 
CD  x 103  C_/Cf 


Number 

ri 

Gertler 

Simple  Formula 

Simple 

4171 

0.00 

0.13 

0.063 

0.17 

4172 

0.05 

0.13 

0.082 

0.22 

4173 

0.15 

0.13 

0.085 

0.22 

4174 

0.20 

0.10 

0.084 

0.22 

‘Based 

on  Rn  ■ 

L 

cf = 

20  x 106  (16) 
2.628  x 10~3 

Computed 

V 


,TSS~~  r 


The  computed  values  of  C as  a function  of  L/D  have  the  same  order 
of  merit  as  the  experimental  values  shown  in  Table  2.  However,  the  computed 
values  of  CR  listed  in  Table*  3,  4,  5,  and  6 and  illustrated  in  Figure  3, 
indicate  that  the  simple  formula  is  relatively  insensitive  to  changes  in  m, 
Cp,  rQ,  and  . 

computed  values  of  C as  a function  of  C , whereas  the  experimental  data 

k p 


For  example,  Table  4 shows  very  little  difference  in  the 

Cp,  whereas  the  experimental  data 
Figure  3 shows  that  except  for  three 


show  a strong  dependence  of  C on  C 

K p 

models  the  experimental  values  of  CR  are  less  than  the  computed  CR  values. 

The  only  series  where  consistent  trends  can  be  seen  in  both  computed 
and  experimental  CR  values  is  for  models  where  L/D  varies  (Table  2) . Both 
experimental  and  computed  CR's  show  decreasing  CR  with  increasing  L/D.  Both 
the  experimental  data  and  the  simple  formula  were  able  to  discriminate  the 
three  worst  bodies  of  this  series;  thus,  it  is  possible  that  the  formula  may 
be  used  to  identify  a really  poor  body  when  L/D  is  being  varied,  but  not 
necessarily  for  variations  of  other  geometric  parameters. 


Models  5242-1,  2,  3 

The  second  series  of  bodies  investigated  was  the  Series  5242  models. 
These  bodies  have  identical  forebodies  with  sterns  of  differing  fullness. 

The  amount  of  parallel  middle  body  on  each  model  was  varied  to  produce  three 
models  of  essentially  constant  volume.  Figure  4 shows  sketches  of  three 
full-scale  hulls  with  changes  in  length  due  to  the  constant-volume 
constraint.  For  the  present  study,  the  hulls  were  treated  as  bodies  of 
revolution  without  appendages. 

Scaled  offsets  were  taken  from  the  drawings  and  polynomials  were  passed 
through  the  offset  points  to  form  a smooth  fit.  The  three  tails  were 
described  by  quartics  of  the  form 

4 3 2 

y = Ax  +Bx  +Cx  +Dx 


The  common  nose  for  the  three  models  was  described  similarly  but  using  a 
quadratic  fit.  Direct  use  of  scaled  offsets  was  also  tried  but  the  rate  of 
change  of  curvature  was  not  smooth  enough  to  avoid  fluctuations  of  the 
pressure  distribution  over  the  body.  Final  pressure  distributions  for  the 
5242  series  can  be  found  in  Figure  8A  of  Appendix  A.  Table  7 shows 
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experimental  results  for  the  Series  5252  models  along  with  the  computed 
values  of  CR/C^  and  the  value  of  CR  calculated  at  a Reynolds  number  of 
20  x 106.  Figure  5 shows  a graphical  comparison  between  the  computed 
and  experimental  values  of  CR. 

TABLE  7 

Residual  Drag  Comparisons  for  the  Series  5242  Models 


Model 

C 

p 

rstern 

Lr/0 

Experimental 
CR  X 103 

Calculated 

CR/Cf 

Simple  Formula 

Calculated 
CR  x 103 

5252-1 

0.674 

2 

0.375 

0.88 

0.23 

5242-2 

0.574 

3 

0.310 

0.076 

0.20 

5242-3 

0.505 

4 

0.255 

0.079 

0.21 

♦Based  on 

Rn^=  20  x 

106  (16) 

The  bare -hull  data  show  that  CR  increases  with  decreasing  LR/D; 

Model  5252-3  with  L /D  = 4 is  the  best  in  the  series.  The  simple  formula, 

K 

however,  predicts  that  Model  5252-2  with  LR/D  - 3 is  the  best  model.  Thus, 

I the  experimental  and  computed  results  are  not  in  agreement  for  these  bodies. 

! Figure  5 shows  that  the  predicted  drag  is  not  very  sensitive  to  change  of 

the  parameter  LR/D.  The  predicted  values  of  CR  are  lower  than  the  measured 

values  of  C , which  was  not  true  for  the  Series  58. 

K 

Kempf  Sterns 

i The  third  series  studied  consisted  of  six  bodies  of  revolution  having 

i hemispherical  noses  and  sinusoidal  sterns  of  varying  fullnesses.  These 

■4  1112 

? sterns  were  investigated  initially  by  G.  Kempf.  ’ The  models  were  4 

j meters  in  length  and  consisted  of  0.3  meter  diameter  pipe  3.1  meters  long, 

with  0.9  meter  afterbodies  described  by  sinusoids  having  different  wave- 

^ lengths.  Sketches  of  the  sterns  are  shown  in  Figure  6. 

a 

i The  point  generation  program  used  the  following  equation  for  each  of 

I the  stern  sections 


19 


r ■ sin  0 ■ sin 


where  0 < 0 < — * 

2 

x is  measured  from  the  aft  perpendicular,  and 
L is  the  length  of  the  rear  section. 

K 

The  tail  radius  of  these  bodies  is  not  zero  at  the  aft  perpendicular,  but 
has  a value  of  37  mm.  To  minimize  possible  and  effects  on  the  pressure 
calculations,  an  imaginary  section  of  0.9  meter  length  and  37  mm  radius 
was  added  to  the  body  points  generated  analytically  (see  Figure  7) . 
Calculations  from  direct  offsets  were  also  tried  but  the  rate  of  change  in 
r was  not  smooth  enough  to  prevent  severe  fluctuations  in  C^. 

Computations  of  the  integral  in  equation  [20]  were  carried  out  only  to 
the  end  of  the  actual  body  and  not  to  the  end  of  the  imaginary  extension. 
Computed  and  measured  drag  values  are  shown  in  Table  8,  a graphical 
comparison  of  CR's  based  on  a model  scale  Reynolds  number  of  20  x 10^  is 
shown  in  Figure  8,  and  final  pressure  distributions  may  be  founded  in 
Figure  9A  of  Appendix  A. 


TABLE  8 

Residual  Drag  Comparisons  for  Kempf  Body  Series 


Model 

lr/d 

Experimental 
CR  x 103 

Computed 

CR/Cf 

Simple  Formula 

Computed 
CR  x 103* 

Simple  Formula 

I 

3.0 

0.07 

0.76 

0.20 

II 

2.5 

0.14 

0.76 

0.20 

III 

2.0 

0.20 

0.78 

0.21 

IV 

1.5 

0.30 

0.081 

0.21 

V 

1.25 

0.42 

0.099 

0.26 

VI 

1.00 

0.75 

0.089 

0.23 

♦Based  on  RnL=  20  x 10  (16) 
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CR  x 10  (EXPERIMENTAL) 

Figure  8 - Comparison  of  Experimental  and  Calculated  CR' 
for  Keapf  Series  Bodies 
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l 
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The  experimental  results  given  by  Kempf  show  a consistent  increase 

in  CR  with  decreasing  LR/D.  With  the  exception  of  Model  VI,  the  simple 

formula  also  shows  this  trend,  although  not  differentiating  nearly  as 

strongly  between  models  as  the  experimental  values.  Model  VI  surely  has 

extensive  separation  which  would  invalidate  the  simple  formula,  since  it 

requires  an  unseparated  boundary  layer.  As  Figure  8 shows,  the  simple 

formula  does  not  adequately  predict  the  very  large  experimental  increase 

of  CD  values  with  decreasing  stern  fineness,  Ld/D. 

K K 

Parallel  Middle  Body  Series 

The  fourth  series  of  body  models  studied  was  based  on  the  experimentally- 

determined  "best"  model  of  Series  58  (Model  4165)  . The  original  equation 

2 t 2 , 3 4 5 6 

y = a, x + a„x  + a„x  + a.x  + arx  + a,x 
1 2 3 4 0b 

with 


a:  = 1.0 

a.,  = 0.837153 

a3  = -8.585996 

a,  = 14.075954 
4 

a5  = -10.542535 
a„  = 3.215422 

D 

was  used  for  the  nose  and  tail  with  increasing  amounts  of  parallel  middle 

body  added  to  generate  a new  series  of  models.  The  experimental  series, 

as  reported  by  Larsen,13  gives  data  for  models  having  30,  40,  50,  and  60 

percent  of  their  lengths  in  parallel -middle-body . 

Sketches  of  the  models  are  shown  in  Figure  9,  computed  pressure 

distributions  are  given  in  Figure  10A  of  Appendix  A,  and  results  of  the  drag 

calculations  are  listed  in  Table  9.  Calculated  values  of  C are  based  on  a 

6 ^ 

model  Reynolds  number  of  20  x 10  , using  Cf  values  calculated  from  the 
16  * 

Schoenherr  formula.  Figure  10  gives  a graphical  comparison  of  the  computed 
and  experimental  values  of  CR. 

16Schoenherr , K.E.,  "Resistance  of  Flat  Surfaces  Moving  Through  a Fluid," 
Transactions,  Society  of  Naval  Architects  and  Marine  Engineers,  Vol . 40 
(1932) . 
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0.2  0.3  0.4  0.S  0.6  0.7 

CR  x 103  (EXPERIMENTAL) 

Figure  10  - Comparison  of  Experimental  and  Calculated  C 's  for  the 
Series  58  Parallel  Middle  Body  Models 


TABLE  9 


■ 


I 


Residual 

Drag  Comparisons 

for  the  Series 

58  Parallel  Middle 

Body  Series 

(Nose  and  tail  shape  fixed,  maximum  diameter  fixed) 

Model 

L/D 

Experimental 
CR  x 103 

Calculated  Calculated 

CR/Cf  CR  x 103 

Simple  Formula  Simple  Formul 

4165 

7.00 

0.07 

0.088 

0.23 

4165-30 

10.00 

0.10 

0.064 

0.17 

4165-40 

11.67 

0.12 

0.055 

0.15 

4165-50 

14.00 

0.14 

0.045 

0.12 

4165-60 

17.50 

0.15 

0.037 

0.10 

Cp  computed  at  RnL  = 20  x 10  . (16) 


The  experimental  data  of  Figure  10  show  that  C values  increase  with 

K 


increasing  amounts  of  parallel  middle  body.  The  simple  formula  shows 


values  of  CR  decreasing  with  increasing  parallel  middle  body.  This  trend 


in  the  computed  drags  is  consistent  with  the  Series  58  data  which  showed  a 
similar  reduction  in  the  values  of  CR  with  increasing  L/D.  It  is  in 
disagreement  with  Larsen's  experimental  data  illustrated  in  Figure  10. 

If  one  considers  a form  of  fixed  nose  and  tail  shape  and  fixed  maximum 
diameter  being  elongated  by  increasing  amounts  of  parallel  middle  body,  it 
would  be  reasonable  to  expect  a thicker  stern  boundary  layer  and  possibly 
larger  values  of  the  residual  (or  form)  drag  due  to  an  increasing  pressure 
defect  at  the  stem.  However,  a comparison  of  the  relative  computed  form 
drags,  given  by  the  product  S-CR,  where  S is  the  wetted  surface  area,  revealed 


the  same  trends  as  given  by  the  values  of  CR  shown  in  Table  9 and  Figure  10. 
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Model  4620-1,  2,  3,  4 Forebody  Series 


The  fifth  series  of  bodies  studied  was  initially  investigated  by 

14 

McCarthy,  Power,  and  Huang,  and  will  be  designated  Series  4620  after 
the  parent  model.  The  parent  model  provided  the  common  tail  form  for  the 
series.  The  series  is  comprised  of  four  bodies  of  revolution  having  bow 
entrance-length/diameter  ratios  (L^/O)  of  0.5,  1.0,  1.82,  and  3.0.  bach 
forebody  has  a prismatic  ratio  Cp^  of  0.667.  The  bluntest  fore-body 


(ll/d 


0.5)  has  a hemispherical  nose.  The  other  three  forebodies  are 


described  by  Granville  polynomials  which  provide  a smooth  junction  with  the 
parallel  middle  bodies.^  These  polynomials  are  defined  by  the  expression: 


y = .833  R(x)  + 10  K^x)  + Q(x) 


where  y 
x 


D 

L 


2V/D 

X/L 

maximum  diameter 
body  length 

dimensional  body  points 


I 


I 

i 

? 

i 


i 


and 

R(x)  = 2 x (x-1)4 
K (x)  = 1/3  x2  (x-1)3 
Q(x)  = 1 - (x-1)4  (4  x +1) 

For  the  series,  the  total  model  hull  volume  was  held  constant  by  varying 
the  length  of  parallel  middle  body.  Sketches  of  the  bodies  are  shown  in 
Figure  11  and  pressure  distributions  for  the  series  are  given  in  Figure  11A 
of  Appendix  A.  The  comparisons  between  the  experimental  data  and  the 
computed  values  of  C^/C^.  are  shown  below  in  Table  10,  and  a graphical 
comparison  of  computed  and  experimental  CD's  is* shown  in  Figure  12. 

Computed  C 's  are  based  on  a model  Reynolds  number  of  20  x 10  . 
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CR  x 10  (EXPERIMENTAL) 

Comparison  of  Experimental  and  Calculated  C 
for  Model  Series  4620-1,  2,  3,  4 


TABLE  10 


Res idual 

Drag  Comparisons  for 

Series  4620-1, 

2,  3,  4 

Model 

Lf:/I) 

Experimental 
CR  x 

Computed 

CR/Cf 

Simple  Formula 

Computed 
CR  x 103 

Simple  Formula 

4620-1 

0.50 

0.25 

0.084 

0.22 

4620-2 

1.00 

0.24 

0.063 

0.17 

4620-3 

1 .82 

0.20 

0.052 

0.14 

4620-4 

3.00 

0.20 

0 . 04  7 

0.12 

‘Based  on 

C^.  computed  at  Rn^  = 20  x 10b. 

(16) 

Both  the  measured  and  predicted  values  of  drag  show  the  .-.are  relative 
order  of  merit  for  this  series.  The  values  of  C increase  with  decreasing 
L^/D.  However,  the  simple  formula  predicts  lower  values  and  a greatei  spread 
between  the  values  of  CR  for  the  four  models. 

OTHER  HULL  MODELS 

A sixth  group  consisting  of  miscellaneous  model  hulls  fitted  with 

polynomials  was  investigated.  Five  models  were  involved:  Model  4620, 

Model  4936,  Model  4627,  Model  5224-1,  and  Model  5290.  Four  of  the  five 

models  were  fitted  with  Granville  polynomials.  Models  4935,  5290  and 

5224-1  were  fitted  using  Granville^'  equations  for  parallel-middie-body 

1 8 

models.  Model  4627  was  fitted  using  Granville  equations  for  forms 

without  parallel  middle  body.  Model  4620  was  fitted  with  a polynomial  of 

Series  58  form,  Y*  = a,  X +■  a X2  + a.  X3  + a,  X4  + ac  X5  ♦ a„  X6  where, 

*2^456 

X = x/L  and  Y = y/D.  The  least  squares  fit  provided  by  these  polynomials 
was  shown  to  be  necessary  to  remove  irregularities  from  drawing  offsets  and 
to  provide  a smooth  distribution  in  C , The  agreement  with  the  offsets  was 
fairly  good.  It  should  be  noted  that  Model  4620  is  an  original  parent  form 
without  parallel  middle  body.  The  4620-1,  2,  3,  4 Model  Series  mentioned 
previously  used  the  tail  for  this  model  but  added  extensive  lengths  of 
parallel  middle  body.  Thus,  computations  for  the  parent  model  are  not 
directly  comparable  to  those  for  the  4620-1,  2,  3,  4 Series. 


Sketches  of  these  model  forms  are  not  included  here,  but  their  pressure 
distributions  are  shown  in  Figure  1 2A  of  Appendix  A.  A comparison  of  the 
predicted  and  measured  drags  for  this  series  is  shown  in  Table  11,  and  a 
graphical  comparison  of  the  CR's  based  on  a Reynolds  number  of  20  x 106 
is  shown  in  Figure  13. 

TABLE  11 

Residual  Drag  Comparisons  for  Miscellaneous  Model  Series 


Experimental 


Computed 

c„/cr 


Computed 
C„  x 103* 


Model 

CR  x 10 

1 

Simple  Formula 

r\ 

Simple  Formula 

4620 

0.10 

0.106 

0.28 

4935 

0.11 

0.083 

0.22 

5 290 

0.15 

0.056 

0.15 

4627 

0.20 

0.089 

0.23 

5224-1 

0.26 

0.061 

0.16 

Based  on 

RnL  = 20  x 106  (16) 

The  experimental  and  computed  values  of  C do  not  agree  as  to  the 
relative  order  of  merit  for  this  group  of  bodies.  Experimentally,  Model 
4620  was  best  and  Model  5224-1  was  the  worst.  The  simple  formula  predicts 
that  Model  5290  is  best  and  that  Model  4620  is  worst.  It  should  be  recalled 
that  the  simple  method  was  unable  to  handle  the  overall  ordering  of  Series  58, 
although  in  that  case,  it  did  correctly  predict  the  worst  body,  which  has 
not  happened  here . 

EXTENDED  STERN  MODELS 

The  last  series  investigated  was  based  on  an  existing  form  represented 
by  Model  4935.  For  this  series,  extended  tail  shapes  were  developed  and 
faired  into  the  original  model.  Model  4935  is  the  original  hull  form.  Model 
4935-2  incorporates  a 10-foot  tail  extension  faired  into  the  original  hull 
approximately  7.5  stations  aft  of  the  parallel  middle  body.  Model  4955-3 
incorporates  an  18-foot  extension  also  faired  into  the  original  hull  about 


«*ra*r*  ’arv mu-  ~,rr: 


f * 


7.5  stations  aft  of  the  parallel  middle  body.  Sketches  of  the  tail  forms 

are  shown  in  Figure  14  and  pressure  distributions  can  be  found  in  Figure  13A 

of  Appendix  A.  Again,  a least-squares  polynomial  fit  was  used  to  fair  the 

model  offsets  in  order  to  obtain  pressure  and  velocity  distributions. 

Previously  unpublished  experimental  values  of  CD  are  shown  in  Table  12 

along  with  the  CR's  computed  from  the  simple  formula.  A graphical  comparison 

of  the  experimental  C values  with  computed  C values  based  on  a model  scale 

Reynolds  number  of  20  x 10  , is  shown  in  Figure  15. 

The  experimental  data  show  a reduction  in  the  values  of  C with 

K 

increasing  tail  length.  The  simple  method  predicts  that  all  hulls  have 
nearly  the  same  value  of  CR.  Thus,  the  simple  method  exhibits  little 
sensitivity  to  changes  in  LR/D  for  these  models.  The  simple  method  was 
also  not  able  to  predict  correctly  the  drags  of  the  Series  5242  bodies 
were  L /D  was  varied,  despite  the  fact  that  the  L /D  changes  for  the 

K R 

Kempf  bodies  were  handled  fairly  well.  The  effect  of  L /D  on  drag  apparently 

R 

cannot  be  satisfactorily  discriminated  by  the  simple  formula  in  all  cases. 


TABLE  12 

kesidual  Drag  Comparisons  for  Model  4935  Extended  Tail  Series 


Stem 

Extension 

Model 

(ft) 

4935-1 

0 

4935-2 

10 

4935-3 

18 

' i 


Experimental 
CR  * 103 

Computed 

CR/Cf 

Simple  Formula 

Computed 

CR  “ 103 
Simple  Formula 

0.23 

0.084 

0.22 

0.17 

0.085 

0.22 

0.13 

0.082 

0.23 

35 


CR  x 10  (EXPERI*NTAL) 

Figure  IS  - Comparison  of  Experimental  and  Computed 
C 's  for  the  Model  4935  Extended  Tail  Series 


rr  1 

i 

CONCLUSIONS 

It  is  apparent  that  the  simplified  method  has  produced  mixed  results. 

The  computed  drags  show  very  little  sensitivity  to  parameters  such  as  C , m, 
r , and  r ^ , as  can  be  seen  from  the  summary  graph  in  Figure  16.  Nor  does  the 
method  appear  to  be  able  to  discriminate  between  bodies  where  many  parameters 
are  varied  at  the  same  time.  In  some  cases  the  method  does  appear  to  be 
capable  of  identifying  order  of  merit  when  parameters  such  as  L^/D  and  L/D 
are  varied.  In  general,  however,  the  simple  formula  cannot  be  recommended 
for  preliminary  design  or  drag  evaluation  purposes.  While  the  formula  does 
in  certain  cases  predict  drag  trends  which  agree  with  experimental  data,  it 
is  by  no  means  capable  of  defining  "optimums'  and  is  very  likely  to  be  in 
i error  for  many  cases.  It  is  difficult  to  make  these  very  broad  conclusions 

more  precise  because  of  the  great  scatter  in  the  results  and  because  of  the  very 
i real  possibility  that  some  of  the  experimental  data  is  incorrect.  More 

detailed  conclusions  require  additional  investigations. 
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APPENDIX  A 


PRESSURE  DISTRIBUTIONS  FOR  ALL  MODEL  FORMS 


Bodies  with  Various  Length/Di ame ter  Ratios,^  (Ref.  15) 
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Figure  5A  - Pressure  Distributions  for  Longitudinal  Flow  for  Series  58 
Bodies  with  Various  Tail  Radii,  r.  (Ref. 15) 


NOTE:  ERRATA  SHEET  TO  REF.  15  STATES  TO  DECREASE  ORDINATE  SCALE  BY  0.05 
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